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Abstract 

• '—I . 

c/3 I Motivated by possible astrophysical and biological applications we calculate 

visible and near UV spectral lines of proflavine (C13H11N3, 3,6-dianiinoacridine) 
in vacuum, as well as its anion, cation, and dication. The pseudopoten- 
tial density functional and time-dependent density functional methods are 

j^ ! used. We find a good agreement in spectral line positions calculated by two 

real-time propagation methods and the Lanczos chain method. Spectra of 
,.^ _ proflavine and its ions show characteristic UV lines which are good candidates 

^ ■ for a detection of these molecules in interstellar space and various biological 



processes. 

Keywords: Astrophysics, Astrochemistry, Biology, Molecular data. 
Density functional theory. Time-dependent density functional theory 



O ' 1. Introduction 



10 ' More than 150 molecules have already been detected in interstellar space 

^ ■ (jj . The number of atoms in by now discovered interstellar molecules is in the 

range from two to more than ten. Broad absorption spectral structures have 
been observed in the near-UV, visible, and near-IR spectra of some stars. 
k> ; They are called diffuse interstellar bands (DIBs) [ity]. Up to now, around 

^ • 300 DIBs have been measured. Although the first such lines were detected 

almost ninety years ago, their carriers are still unknown. Many DIBs carriers 
have been suggested, but two main groups are processes on the astrophysical 
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dust grains and molecular transitions. The molecular origin of DIBs is more 
dominant hypothesis in the recent literature. Various simple or complex car- 
bon containing molecules have been suggested as possible DIBs carriers. For 
example, small linear carbon chains, fullerenes, and especially polycyclic aro- 
matic hydrocarbons (PAHs) have been investigated as candidates for DIBs 
carriers j#^8|. Although various neutral and ionized PAHs, small and large 
ones, have been extensively studied, the analysis still did not unambiguously 
confirm agreement between any specific PAHs spectra and DIBs. A few nar- 



row DIBs have been observed in the Red Rectangle nebula [91 



rr^i 
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Substructures 
Features in 



in prominent bands of some narrow DIBs have been seen [id . 
these resolved profiles are consistent with rovibrational transitions in organic 
molecules. It is also important to study various molecules where atoms of 
other chemical elements substitute C and H atoms in PAHs. Nitrogen is a 
very important component of biological materials, and it is one of the most 
abundant chemical elements in the Universe. Several molecules that contain 
nitrogen have already been detected in meteorites, interstellar dust particles, 
and Space, or they have been analyzed by computational methods 12Nl4 . 



Proflavine, CiqH 



N3 (3,6-diaminoacridine), is a molecule well known in 



molecular biology |l5Hl7l|. This molecule interacts with nucleic acids, inter- 
calates into DNA, and blocks activities of many viruses. C13H11N3 is a typical 
example of amino- acridines which have been tested and used as anticancer 
drugs. It has been suggested that proflavine is one of molecules which may 



act as "midwifes" [17| . It is possible that this type of molecules increased a 
probability of DNA and RNA synthesis when biological matter was formed. 
For example, proflavine increases thousand times a kinetics of the synthesis 
of two small DNA or RNA molecules into larger one. Therefore, proflavine 
is one of the candidates for molecules which induced the origin of life on 
Earth 18|. Optical spectra present an efficient way to detect proflavine, its 
ions, and derived nanostructures in biological materials and processes. In 
a comparison with PAHs (for example oligoacenes) , available experimental 
and theoretical data for proflavine are scarce. In ref. 3 it has been suggested 
that dyes are good candidates for studies of a DIBs problem and, as a sup- 
port of this idea, it was stated that proflavine absorbs at 445 nm in water 
at pH 7. However, existing measurements of proflavine spectra are done in 
aqueous and other solutions [l6|, Il9l-l22|. It has been found that the optical 



spectrum of proflavine is very sensitive to the type of a solvent and pH factor 



19l |2l| . There is a possibility that spectra of proflavine and its ions in vac- 



uum and under astrophysical conditions are substantially different from ones 



in solutions. In addition, it is known that many materials from the group of 
organic dyes show strong lines in a visible spectrum. For example, absorption 
lines of basic dyes such as methylene blue, thionine, neutral red, safranine, 
acridine orange, acridine yellow, and proflavine are all in an interval from 
400 A to 800 A 20|. Although many organic dyes are good candidates for 



DIBs studies, because of a "midwife" idea [17[ and possible applications in 
astrobiology, we have selected proflavine for our study. 

It is known that many molecules in interstellar space exist in its ionized 
state. It has been found that positive and negative ions of PAH molecules 



play an important role in the chemistry of diffuse interstellar clouds [23 
Therefore, it is also important to analyze ions of proflavine. We use pseudopo- 
tential density functional (DFT) methods to compute the ground-state prop- 
erties of proflavine and its ions, as well as the pseudopotential time-dependent 
density functional theory (TDDFT) to calculate their optical spectra. Pseu- 
dopotential methods are computationally less demanding than all-electrons 
methods and offer a possibility to study larger systems. Methods based on 
the pseudopotential DFT are now substantially developed to provide rehable 
results for structural, electronic and optical properties of various materials 



24( 1 . TDDFT |25| is recently used to describe optical properties of many 
materials p-lsl, |26|-|29||. It is known that electronic spectra of PAHs are in 
the UV and visible part of the spectrum. Optical spectra of PAHs and many 
other organic molecules in this region of the spectrum are well described by 



TDDFT tern, 27-29 



2. Computational methods 

The ground-states of all proflavine related molecules are optimized in- 
dependently using the Quantum ESPRESSO program package [SO"]. For a 
majority of optical spectra we use a real-time pro pag ation method within 



TDDFT as implemented in the Octopus code |26l . |31[ . The first step is to 
calculate the ground state. Then, the system is perturbed by a small time- 
dependent potential SVextif, t) = —Kz5t. This perturbation is turned off, and 
the time-evolution of the system is followed. If the dipole momentum of the 
system is d{t) = —e J p{f, t)fdf, then its dynamical polarizability is obtained 
from a{u) = - J e*'^*[(i(t) — d{0)]dt. Optical spectrum is calculated from 

S{u) = —Ima{uj). (1) 

TT 



The real space method and the spacing of 0.13 A are used. The simulation 
box is formed by adding spheres of the 4 A radii around each atom. The 
time-dependent equations are integrated using the time step of 6t = 0.0012 
h/eV. The time evolution of ~ 7.9 fs is followed. The approximated enforced 



time-reversal symmetry algorithm is applied [32 



We use the TDLDA approximation and the Perdew-Zunger exchange- 
correlation functional [33| because of its good numerical stability in the pseu- 
dopotential TDDFT method. The Troullier-Martins pseudopotentials are 



used J3J]. It is known that hybrid functionals in all-electron DFT methods 



produce results which are often in better agreement with experiments than 



those obtained using LDA and GGA [35|. In contrast to all-electron codes, 



applications of hybrid functionals in the plane wave pseudopotential DFT 
methods are still in development. This is even more pronounced for pseu- 
dopotential TDDFT methods. We compare structural parameters of neutral 
proflavine calculated in LDA with those obtained using the BLYP functional 
in the Quantum ESPRESSO |36l-l38| . We also calculate optical properties 



of neutral proflavine using the B3LYP functional in the latest version of the 
Octopus code, where hybrid functionals are still in development, and find an 
agreement between LDA and B3LYP optical spectra. However, because of 
much better numerical stability, as well as availability and reliability of corre- 
sponding pseudopotentials, we present results for the TDLDA functional. It 
has been shown that the TDLDA pseudopotential approximation in a time- 
dependent evolution gives optical properties in agreement with experiments 
and other calculations, even for complex biological molecules such as the 
green fluorescent protein 27| and DNA bases |28]. 

We also compare optical spectrum of proflavine calculated by the Oc- 
topus code with one obtained by the TDDFT module |29|, |39| of Quantum 
ESPRESSO J30|- This module presents a new algorithm for the time depen- 
dent density functional perturbation theory. A susceptibility is determined as 
an off-diagonal matrix element of the resolvent of the Liouvillian superopera- 
tor. Computations are done using a Lanczos continued fraction method. The 
optical spectrum of proflavine is calculated using the ground state computed 
in the Quantum ESPRESSO package. Periodic boundary conditions and the 
tetragonal simulation cell of sides length 30 a.u. and 16.7 a.u. are used. A 
local density Perdew-Zunger exchange-correlation functional |33|, and ultra- 
soft pseudopotentials of the RRKJ type [40| are applied. The kinetic energy 
cutoff of 25 Ry for wavefunctions, and 200 Ry for the charge density are 
used. The convergence properties in the Lanczos method are checked and 




Figure 1: The structure and electron density in the plane of a neutral proflavine molecule. 
Three nitrogen atoms are labeled by "N". Large balls represent carbon atoms, whereas 
small ones are hydrogen atoms. Isosurfaces are shown for 0.3 e/ajj isovalues, where qq is 
the Bohr radius, and e is the electron charge. The figure is prepared using the XCrySDen 
package [41 1. 



the optical spectrum for 3500 recursion steps, and with the broadening of 
0.01 Ry is shown in this work. 



3. Results and Discussion 

The structure and electron density of proflavine molecule is shown in 
Fig. [H It is found that different exchange functionals and corresponding 
pseudopotentials produce results for geometry parameters which are suffi- 
ciently close to each other. For example, in proflavine molecule average 
N-C distances are 1.365 A for LDA and norm-conserving pseudopotentials 
0, 1.349 A for LDA and ultra-soft pseudopotentials jio], 1.371 A for the 



BLYP functional |36l-l38| . Average C-C distances are: 1.395 A for LDA and 
norm-conserving pseudopotentials, 1.396 A for LDA and ultra-soft pseudopo- 
tentials, 1.413 A for the BLYP functional. This should be compared with 
computationally generated data available in the NIST database [42]: 1.365 
A for average N-C, and 1.406 A for average C-C distances. 

It is calculated that the LDA Kohn-Sham HOMO-LUMO gap for neutral 
proflavine is 2.240 eV. The quasiparticle corrected HOMO-LUMO gap is 
obtained using two formulae JTl, |43|, |4^. From QPgap = -^'a^+i + E^^i — 2E]\f, 
where E^ is the total energy of N-electron system, the value of 4.509 eV is 
obtained. The second formula for the quasiparticle corrected HOMO-LUMO 



gap QP^ 



2 
gap 



-N+1 



- f^ 



where e] is the i-th eigenvalue of the j-electron 



system, produces the value of 4.511 eV. 



Optical spectra of proflavine and its ions are shown in Fig. [21 It is found 
that the overall spectra are rather similar to each other. However, differences 
in the position and shape of lines exist. Optical transitions in ions occur at 
lower energies than in a neutral proflavine molecule. The same relation be- 



tween spectra has been found for neutral and charged PAH molecules |45 . 
We find that the first pronounced absorption peak in spectra of proflavine 
and its ions appears at ~ 3 eV. The first absorption peak in PAH anthracene 
(with the same number of cycles) is at 5 eV (7| . Figure |3] compares spectra 
of neutral proflavine calculated by two exchange-correlation functionals in 



the Octopus code [3l| with the same result obtained by the TDDFT mod- 



ule J29|, |39| of Quantum Espresso [30j • Although numerical algorithms and 



pseudopotentials used in these methods are different, the agreement between 
shapes of the spectra and positions of lines is good. It has been measured 
that a position of a visible absorption line of proflavine in water at room 
temperature changes from 444 nm to 394 nm when pH changes between 7.0 



and 14.0 21 . These authors concluded that the value of 394 nm corre- 



sponds to the neutral form of proflavine in water. The same line is at 397 



nm for proflavine in purified toluene [2l[. The value for the maximum of the 



absorption spectrum of proflavine in water reported in ref. [22[ is 400 nm. 

The comparison of calculated lines and corresponding DIB lines 46|-|48 
is presented in Table [H Data for a visible spectrum show that only line 
of a neutral molecule is within FWHM of known DIB lines. Positions of 
the UV line nearest to the visible spectrum, as well as of a characteristic 
strong line which appears in all spectra between 8.35 eV and 8.55 eV, are 
also shown in Table [H Many organic molecules exhibit absorption lines in 
the UV spectrum. Therefore, the measurements and identification of UV 
spectral features are important for organic materials in Space. For example, 
the interval between 1320 A and 1440 A in the UV spectrum measure d by 



the Space Telescope Imaging Spectrograph (STIS) has been analyzed J49 
The Cosmic Origins Spectrograph (COS), recently installed on the Hubble 
Space Telescope, will provide new data for spectra of organic molecules in 
UV [4]. Calculations based on a model of the state of hydrogenation and 
the degree of ionization in the diffuse interstellar medium have shown that 



small PAHs with 15-20 carbon atoms are destroyed in the UV field [50 
Although proflavine and its ions are rather small, there is a possibility that 
they exist in Space adsorbed on suitable grains, or as intercalators in some 
bigger molecules. That could even change their optical spectra and give line 
positions in better agreement with DIBs. 
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Figure 2: Optical absorption spectra: (a) C13H11N3, (b) CiaHnNa^i, (c) CisHnNa+i, 
(d) Ci3HnN3+2. 
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Figure 3: Spectral lines of neutral proflavine in the near UV (Ai) and visible (A2) re- 
gion calculated by three pseudopotential TDDFT methods: (a) the TDLDA real-time 
propagation method JsJ, Ai = 277.37 nm, A2 = 407.84 nm, (b) the TDB3LYP real-time 
propagation method [31|, Ai = 280.54 nm, A2 — 424.66 nm, (c) the Lanczos chains method 
iillfllli, Ai = 271.76 nm, A2 = 412.76 nm. 



Table 1: Comparison of calculated optical spectral lines of proflavine and its ions with 
the closest known DIB lines 461448 1. FWHM of measured lines is also shown, as well 



as positions of two UV lines to facilitate a comparison with future measurements of the 
ultraviolet range in spectra of interstellar organics. 

"Molecule Lines (A) DIB (A) FWHM (A)~ 

C13H11N3 4078.42 4066.0 [46] Tsl) 

2773.69 
1450.11 
C13H11N3-1 4305.00 4363.86 [471 0.46 



2883.35 

1484.84 
Ci3HnN3+i 4217.15 4175.46 [48] 17.2 

2805.07 

1455.22 
C13H11N3+2 4232.08 4259.01 [47] 1.05 

2831.05 

1458.64 



4. Conclusions 

We calculated optical spectra of proflavine, its anion, cation, and dication 
using pseudopotential density functional theory methods. A comparison of 
three TDDFT pseudopotential methods shows that they produce spectra of 
proflavine which qualitatively agree well, and line positions differ 0.14 eV, at 
the most. The difference of an average result (calculated by three TDDFT 
methods) for a visible line position and a corresponding average experimental 



value for a neutral proflavine molecule in water [21|, |22| is also 0.14 eV. The 
shapes of spectra for proflavine and its ions are rather similar to each other, 
but the change in the line positions induced by a different charge is apparent. 
UV lines are characteristic for spectra of proflavine and its ions and give an 
efficient way for their detection, either in Space, meteorites, astrophysical 
dust particles, or in biological processes. However, the UV spectral region is 
much less studied than visible one. It is important to analyze available data 
and to perform new studies of substituted polycyclic aromatic hydrocarbons. 
We hope that our work will generate an interest for measurements of optical 
spectra of proflavine and other organic dyes in the gaseous phase and under 
various astrophysical conditions. 
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